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We consider the role of nonlamellar-forming lipids in biological membranes by examining fluctua- 
tions, within the random phase approximation, of a model mixture of two lipids, one of which forms 
' lamellar phases while the other forms inverted hexagonal phases. To determine the extent to which 

^s^j , nonlamellar-forming lipids facilitiate the formation of nonlamellar structures in lipid mixtures, we 

examine the fluctuation modes and various correlation functions in the lamellar phase of the mix- 
ture. To highlight the role fluctuations can play, we focus on the lamellar phase near its limit of 
stability. Our results indicate that in the initial stages of the transition, undulations appear in the 
lamellae occupied by the tails, and that the nonlamellar-forming lipid dominates these undulations. 
The lamellae occupied by the head groups pinch off to make the tubes of the hexagonal phase. 
Examination of different correlations and susceptibilities makes quantitative the dominant role of 
the nonlamellar-forming lipids. 
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I. INTRODUCTION 



It is well known that the biological, lipid bilayer membrane contains a mixture of lipids some of which, alone in 
aqueous solution, do not form bilayers at alL This presents the interesting question as to what possible function in 
q ■ the membrane these lipids can be servingoEl. There are two lines of thought on the answer. The first is that the 
'— nonlamellar-forming lipids, roughly characterized as having a small head group and long, splayed tails, can with these 
tails increase significantly the local lateral pressure distribution in the interior of the bilayer. This increased pressure 
may be necessary for a protein or peptide to functioncTEl. The second is that the nonlamellar-forming lipids can 
I . make less expensiije the energetic cost to form a small, nonlamellar region which may be needed in processes such as 
\Q ' membrane fusionLrU. Experiment shows that there is the expected correlation between the presence of nonlamellpi- 
CNJ . forming lipids and the ability of the lamellar phase to make a transition to a nonlamellar, inverted-hexagonal oneuE3. 
There is also a correlation betweeiythe concentration of such lipids and the ability to make other nonlamellar structures, 
such as occur in membrane fusiorJij. ■— ■ 

In an attempt to shed some light on the role in membranes of nonlamellar-forming lipids, we first introducedlia a 
relatively simple and tractable model of a lipid. We showed that it exhibited the polymorphism observed in experiment 
as a function of a single architectural parameter, the fraction of the total lipid volume occupied by the head group. 
We also showed that the phase diagram of a model aqueous lipid whose head group fr actio n was similar to that 
of dioleoylphosphatidylethanolamine was in good agreement with the measured diagramPH 1 !. We then considered 
a mixture of two lipids with the same headgroup but different tail lengthE£l. Thus they were distinguished by the 
different relative volume fractions of their headgroups such that one formed lamellar, L a , phases while the other 
formed inverted-hexagonal, Hu, ones. We examined the density profiles of the tail segments from the two different 
lipids in each of the two phases. In particular we determined that, in the Hu phase of the mixture, the density of 
the tail segments of the nonlamellar-forming lipid varied around-the Brillouin zone, with its density being largest in 
the direction between next-nearest-neighbor cores, as expectedtl We were able to make quantitative this expected 
variation, and found the rdLative^difference in tail segment densities to be on the order of a few percent. 

The calculations in Ref£j andlij were carried out within self-consistent field theory, (SCFT), and therefore ignored 
fluctuations. In this paper, we consider the effect of Gaussian fluctuations about the SCFT solution. This enables us 
to determine the effect of fluctuations, within Gaussian order, on the transition from the L a to the Hu phase, and to 
observe the initial stages of the path between them. In order to highlight the effect of fluctuations, we examine the 
metastable lamellar phase near its spinodal in a region of the phase diagram in which the Hu phase is the stable one. 
We find that undulations occur in the lamellae occupied by the tails, and that the lamellae occupied byJdie head groups 
pinch off to form the tubes of the Hu phase. This is similar to scenarios proposed earlier by Hui et al.EZI and CaffreyEj 
and which arose from their experimental observations. The hexagonal phase that is initiated by the lowest energy 
fluctuation mode is characterized by a lattice parameter which is within 1% of the stable Hu phase, and has the same 
orientation: one in which the tubes are coplanar with the lamellae from which they were formed. This coplanarity 
is in agreement with experimentlla. We observe that the nonlamellar-forming lipids dominate the undulations in the 
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tail region, while the lamellar-forming lipids tend to fill in the regions between the cylinders as they pinch off. To 
make quantitative the effect of the nonlamellar-forming lipids, we examine three different measures of their role in the 
transition; the correlation between the density fluctuation and the order parameter, the ratio of the susceptibilities of 
the two different lipids to fluctuations at the critical wavevector, and the overlap of the order parameter fluctuations 
of each lipid with the equilibrium Hu structure. All these measures show the nonlamellar-forming lipids to play a 
significant role in bringing about the transition. 

Having calculated the fluctuations in Gaussian approximation, we are able to calculate all structure factors. We 
show here the order parameter-order parameter structure factor, which is the one most readily measured, as well as 
the density-density structure factor of either lipid. r—. 

In the next section we review the theory of Gaussian fluctuations about an ordered phaseE]~E3 and extend it to our 
model of lipid mixtures. Results are presented next, and we conclude with a brief summary. 



II. THEORY 



We consider an anhydrous mixture of n% lipids of type 1 and lipids of type 2 in a volume V. Below we shall 
choose their architecture so that type 1 lipids form lamellar phases while type 2 lipids form Hu phases. All lipids 
consist of the same head group of volume Vh and two equal-length tails. Thus we model mixtures of lipids drawn 
from a homologous series, such as the phosphatidylethanolamines studied by Seddon et alE3. The local density of 
the headgroups of lipid L = 1,2, measured in units of the density 1/vh, is denoted <J>/ l; i(r). Each tail of lipid L 
consists of NctL segments of volume Vt- The local density of lipid tails of lipid L, again measured in units of l/vh, is 
denoted "^^(r). The local volume fraction of these tail segments is (2Nvt/vh)$t,L = lt^t,L- The sole architectural 
parameter, /l = Vh/(vh + 2NaLVt) = 1/(1 + ctLjt), which characterizes each lipid is the relative volume fraction of 
its headgroup. 

The only interaction in the system is between the headgroups and the tails. The interaction energy E is 

1 2N\ f 
— E[$ hA + $ t)1 + $ t , 2 ] = / [* M (r) + $ h ,2(r)][* M (r) + $ t ,2(r)]dr, (1) 

Vh J 

where \ ' IS the strength of the interaction, and T is the temperature. 



The partition function of the system can, without approximation, be writtenli3ui 

Z=f[ [ P$ ft , i W ft , £ OTt )£ OW t , £ e X p[-n/fcT], (2) 

T 1 ** 



where T> denotes a functional integral, and the grand potential VL is given by 

fi = -^-^2^ L Q L [W h , L ,Wt,L] + J dr[W h ,L$h,L + Wt,L$t,L}^ 

+ J5[$ M + * fc|3> 4t 1 i + *t,2], (3) 

with Q.L[Wh.L, Wt,L] the partition function of a single lipid of type L in the external fields Wu.l an d Wt l- 

The integrals over the densities &h,l and $>t,L could be carried out as the densities appear only quadratically, but 
the integrals over the fields Wh,L and Wt,L can not be. Thus some approximate evaluation must be made. The SCFT 
consists in replacing the exact free energy, —kTlnZ by the extremum of Q. In addition, we cxtrcmizc this potential 
subject to an incompressibility constraint that the sum of all head and tail densities is unity everywhere. We denote 
the values of the §h,Li Wh,L: ®t,L an d Wt,L which extremize f2 subject to this constraint by lower case letters. They 
are obtained from the following set of self-consistent equations: 

SQr 

<t>h, L {r) =-z L x , , £=1,2, (4) 



'Sw h<L (r) ' 

SQl 
' 5w t ^ L (r) 1 



',i( r ) = ^7%T. L = l,2, (5) 



«flh,i(r)=2xiV£;&,j/(r)-£(r), £ = 1,2, (6) 
v 

w ttL (r)=2 X Nj2<t>h < L'(r)-ltZ(r), L = 1,2, (7) 
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l = ^<M,L(r)+7tX^( r )> 

L L 



(8) 



where £(r) is the Lagrange multiplier used to enforce the incompressibility constraint. The free energy in this 
approximation, fl sc f, is obtained from Eq. ^ by replacing the densities and fields by their self-consistent values. 

To include fluctuations about this self-consistent field result, one decomposes the fields, Wh,L and Wt,L, and the 
densities, §h,L and &t,L, into their self-consistent field values plus a fluctuating part; Wh,l = Wh,L + 8Wh,L, etc. One 
then substitutes these decompositions into the free energy, Eq. ||. Terms linear in the deviations vanish by definition 
of the self-consistent field approximation. We consider all quadratic terms in this paper so that the fluctuations are 
accounted for within Gaussian approximation. Thus the exact partition function of Eq. is approximated by 
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Z » exp[-n scf /kT] J] / V5$ h , L V5W h , L -D5$ t , L V5Wt,L exp[-n^/kT}. 

L=l J 



(9) 



The expression for ftP) is most easily written in matrix form. Introduce the column matrices 



sw = 



( sw hA 

SWt,i 
5W h , 2 

\ sw< 2 



(10) 



and 



5^> 




(11) 



and the square matrices 



J = 



10 1 



10 1 





(12) 



K 



( ziChh.i 


V o 



z\C t t,i 





ZlChh,2 Z 2 Cht,2 

z%Cth,i z-iCtt,i 



where 



C a p, L {r,r) = V— j-^j^ '— a,/3 = hort. 



8W a . L {v)5W fj , L {v') 

In terms of these matrices, the second order correction, , to the self-consistent free energy is 



VkT V 



drS® 1 (r) J5$(r) 



- i J drW T (r)<5$(r) - ^ J drdr' 8W T {v)K{v, r')SW(r'), 



(13) 



(14) 



(15) 



where the superscript T denotes transpose. The Gaussian integrals in Eq. || over the four fields SWh,L, SWt,L, 
L = 1,2, can now be carried out. This will leave W 2 ' a functional of the quadratic density fluctuations. These 
fluctuations are not integrated over in order to leave displayed their coefficients, which are the inverse of the desired 
density-density correlation functions. It is the inclusion of the field fluctuations to Gaussian order for given density 
fluctuations which constitutes the random phase approximation. 
Integration over the field variables in Eq. ^| yields 



3 



Z«M>exp[-Q sc/ /fcT] Y[ [ V5$ hiL V5$ t!L exp[-n rpa /kT}, 

L=l"' 

where the factor J\f is the reciprocal of the square root of a determinant which is of no interest here, and 

kT 



(16) 



n rpa = ^—^ J dvdv'5® 1 (r) [A X NVJ5{r - r') + K (r, r')J <f*(r / ). (17) 

In addition to the random phase approximation, we now impose the constraint of incompressibility, which reduces the 
number of fluctuating fields from four, 5&h,L, S®t,L, L = 1,2, to three. We choose these three to be natural order 
parameters. Define the total local volume fractions of each lipid 

^L = ^h,L+lt^t,L, £ = 1,2, (18) 
and the local difference between head and tail densities of each lipid 

*L = $h,L-— $t,L, L=l,2. (19) 

Note that as defined, the integral of ^f L over the whole system vanishes. In terms of these densities 

$h,L = h$L + (1 - (20) 
$t,L = a L f L [<f> L -* L ]. (21) 

We choose the three independent fluctuating quantities to be S^>i and 5 1 $>2, the fluctuations in the difference in head 
and tail densities of each lipid, and (5$i, the fluctuation of the total density of lipid 1. Because of the incompressibility 
constraint, <5$i = — <5<j> 2 . The reduction in variables is easily written in terms of the column matrix 



(56 = 5*! (22) 



and the 4x3 matrix 





U= - u , (23) 



so that 

(5$ = USB. (24) 
Then the correction, tt rpa to the self-consistent free energy can be written 

kT f 

n rpa = drdv'Se T (r)U T [A X NJ5{r- r') + ^- 1 (r,r')] USQ(r'), 

■ J drdr'Se T (r)K~ p 1 a (r 7 r')Se(r / ). (25) 



kT 
2v h V 



There remains the calculation of the partition function Q,L[Wh,L,W t: L] of a single lipid in the external fields Wh,L 
and Wt.L- Before carrying this out, we must specify further the way in which the lipid tails are modeled. They are 
treated as being completely flexible, with radii of gyration R g x = (iV^a 2 /6) 1 / 2 for each tail. The statistical segment 
length is a. The configuration of the i'th lipid of type L is described by a space curve rj^s) where s ranges from 
at one end of one tail, through s = at which the head is located, to s = aj,, the end of the other tail. 

Because the tails are completely flexible, one can define the propagator gi,(r, s\r') which gives the relative probability 
of finding the zero'th segment of a tail of a type L lipid at r' and the segment s of the same tail at r: 



/Ti,i,(s)=r r r-s 

q L {r,s\r') = / Vr ltL (s)exp dt W t . L (y LL (t)) 

Jr UL (o)=r' L Jo 



0<s<ai/2. (26) 
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In this expression, Vy^l(s) denotes a functional integral over the possible configurations of the tail of the lipid of 
type L and in which, in addition to the Boltzmann weight, the path is weighted by the factor V[r^L] 0, s] with 



V[r, 81,82] = A/"exp 



ds 



dr(s) 2 



SB*J B1 1 ds 



(27) 



where Af is an unimportant normalization constant and R g = (Na 2 /6) 1 / 2 is the radius of gyration of a tail of length 
N. This propagator is of central importance, both within the SCFT and the random phase approximation. Because 
the tails are flexible and execute a random walk in the presence of the potential Wt^L, the propagator can be obtained 
from the solution of the diffusion equation 

dqL(r ds SW) = 2^V V(r, 8|r') - W t , L (r)q L (r, s\r'), (28) 

subject to the initial condition 

<?i(r,0|r') =<5(r-r'). (29) 
The integrated propagator, or endpoint distribution function, 

q L (r,s) = J dr'q L (r,s\r') = J dr'q L (r', s|r), (30) 

is also useful. 

From the propagator, the partition function of a single lipid of type L is obtained via 

QL[W h , L ,Wt, L } = J dr 1 dr 2 dr 3 q L (r 3 ,a L /2\r 2 )exp[~W h .L(r2)}qL(r2,a L /2\r 1 ). (31) 
Functional derivatives, such as those required in Eqs. |J, |^, and [l4|, are obtained from the above and 

AM/ r f r ? =~ f dtq^s-AvMr^). (32) 
The second functional derivatives of the single lipid partition function are expressed in terms of the propagator; 

C hh , L (r,r') = VS(r~r')q L (r,a L /2)e- w '^q L (r,a L /2) : (33) 

C ht , L (r, r') = 2Ve- w '^q L (r, a L /2) f ' ds q L (r, a L /2 - s\r')q L (r', s), (34) 

Jo 

C thjL {v\v) = C ht<L {v,v'), (35) 

Cu,lM = 2[F L (r,r')+F L (r',r) + G L (r,r% (36) 



where 



and 



rai/2 

F L (r,r') = V ds dt q L (r, a L /2 - s)q L (r, s - t\r') x 

Jo Jo 

J dr 2 q L (r', t\r 2 )e' w ^ r ^q L (r 2 ,a L /2), (37) 

G L (r,r') = V ds dt q L (r, a L /2 - s)q L (r', a L /2 - t) x 

Jo Jo 

dv 2 q L (r,s\r 2 )e- w ^^q L (r 2 ,t\r'). (38) 



Although Eqs. |26j-|38j are exact for arbitrary fields Wh^i?) arid Wt,L( r ), only the mean- field solutions Wh,L{^) and 
u) t> i(r) are needed in the random phase approximation. 

Correlation functions in real space can now be calculated. Because there are three independent order parameters, 
defined in Eqs. |l^ and [ll], there are six independent correlation functions, < 5<di(r)5@j (r') >, i, j — 1,2,3, where the 
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brackets indicate an average using the partition function evaluated within the random phase approximation. We will 
be interested in correlation functions which are various linear combinations of these, so we define a general fluctuation 

S9 a (r) ee oi5$i(r) + a 2 5*i(r) + a 3 5V 2 (r), 
= (01 a 2 o 3 ) <5*i (r) 

= a-Se(r), (39) 

and consider the correlation of two such general fluctuations a T - < SQ(r)5Q T (r 1 ) > b. When the ensemble average 
is carried out, one obtains 

a T - < Se(r)Se T (r') > b = V ± a T • K rpa {v, r') • b, (40) 

where K rpa (r, r') is the inverse of K~p a (r,r') defined in Eq. ^5| 

Because the phases of interest to us are ordered, it is advantageous to exploit the space-group symmetry and, 
following Shi et alB3, to expand all functions of position in terms of Bloch states; 

VWr) = e lk - r ^ U „ k (G) e jG - r . (41) 

G 

Here k is a wavevector in the first Brillouin zone, n is the band index, the vectors G are the reciprocal lattice vectors 
of the ordered phase, and the wavefunctions u„k(G) satisfy the Schrodinger equation 

2R 2 g (k + G) V lk (G) + W t {G - G>„ k (G') = e„ k u„ k (G). (42) 

G' 

The potentials, W t (G), are just the coefficients of the expansion of the periodic field w t (r), Eq. [7] 

w t ,i(r) =w t , 2 (r) =^e 4G ' r ^(G). (43) 

G 

The expansion of functions of two spatial coordinates K(r,r'), such as the correlation functions, takes the form 

K(r,r')= £ W „ „ ^ k (r)C'k'(r') (44) 

nk;n'k' 

Experimentally measured structure factors are Fourier transforms of the correlation functions 
5 ah (q) ee a T ■ i J dvdv'e~^- r '^ < <56(r) ( 5e T (r / ) > -b 

= e t ,^J2\K rpa \ ^ , WG (GK, q _ G (G)-b, (45) 

V — * L J nq-G:n'q-G 



n.n' 



where q — G lies in the first Brillouin zone. 

All calculations are carried out in the basis of Bloch functions. Therefore the important propagators (?i(r, s\r') 
must be expanded in them, and the procedure outlined above in real-space is followed in the, space of Bloch functions. 
This is straightforward, but tedious. Many of the details can be found in the Appendix o£3, and we shall not repeat 
them here. We turn, instead, to our results. 



III. RESULTS 



We have chosen to model two lipids with the same head group, but with different length tails; lipid 1 characterized 
by ax = 1 so that its tails are of length N, and lipid 2 characterized by a 2 = 1.5 so that its tails are of length 1.5 
N. Because the headgroups are identical, the ratio 7 t ee 2Nv t /vh is the same for each lipid, and we have chosen 
7t = 2.5. With these parameters, the volume of the head group relative to that of the entire lipid is, for lipid 1, 
fi ee 1/(1 + ai7t) = 0.2857, while that of lipid 2 is f 2 = 1/(1 + a 2 j t ) = 0.2105. For comparison, the relative head 
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group volume of dioleoviBhosphatidylethanolamine calculated from volumes given in the literature^ is / = 0.254. 
From our previous worked, we know that anhydrous lipid 1 forms a lamellar phase, while anhydrous lipid 2 forms an 
inverted hexagonal phase. 

For [Orientation, we reproduce in Fig. 1 the phase diagram of the anhydrous system of the two lipids calculated 
earlienl3. The temperature, T* , is defined in terms of the interaction strength T* = \/2xN\, and 8 is the volume frac- 
tion of lipid 1, the lamellar- forming lipid. Although we know of no phase diagrams for anhydrous mixtures of lamcllar- 
and non lamellar- forming lipids with the same headgroup but different length tails, the system for which we have 
made our calculation, there are results for mixtures of lamellar-forming phosphatidylcholine and hexagonal-forming 
phosphatidylethanolamine. Our phase diagram is similar to those observed in anhydrous mixtures of dilinoleoylphos- 
phatidylethanolamine and palmitoyloleoylphosphatidylcholiptpj, and for mixtures of dioleoylphosphatidylcholine and 
dioleoylphosphatidylethanolamine and 10% water by weightEH. Each show a significant region of cubic phase between 
the inverted hexagonal phase, which dominates at low concentrations of the lamellar-forming lipid, and the lamellar 
phase, which dominates at high concentration. The dashed line in Fig. 1 is the calculated spinodal of the lamellar 
phase; that is, for volume fractions of the lamellar-forming lipid which are to the left of (i.e. less than) this line, the 
lamellar phase is absolutely unstable. To the right of this line, but in the regions in which either the inverted hexago- 
nal, Hn or inverted gyroid, Gu, is stable, the lamellar, L ai phase is metastable. At high temperatures, the system is 
disordered, D. We have indicated on this phase diagram the two points at which we have calculated structure factors; 
very near the lamellar spinodal line at T* — 0.0625, 8 = 0.356, at which point the fluctuations will be large, and at 
T* = 0.0500, 8 = 0.469, far from the spinodal, where fluctuations will be much reduced. At the former point, the 
stable phase is Hu, while at the latter it is the Gu. 

The structure factor most readily measured experimentally is the Fourier transform of the order parameter-order 
parameter correlation function < 5 , J , (r)<5 , i , (r / ) >, with <$\I>(r) = <5^' 1 (r) + S"i>2( r ) the total order parameter. This order 
parameter is proportional to the difference in the densities of headgroups and tails. The structure factor is obtained 
by setting the vectors a = b, defined in Eq. to be (0, 1,1). It has been calculated at T* — 0.0625, 8 = 0.356, and is 
shown in Fig. 2. The normal to the lamellar planes is in the z direction. The wave vectors k z and k x are measurecLin 
units of 2tt/D, where D is the equilibrium lattice spacing of the lamellar phase. As in the block copolymer systemEil, 
there is a large response about the positions of the Bragg peaks (k x ,k z ) — (0,±1), (0, ±2), and (0,±3). In addition 
there are four satellite peaks at (±0.846, ±1/2). As shown below, these fluctuation modes lead to undulations in the 
planes of head-groups causing them to pinch off, and to form tubes. The array of tubes produced by these fluctuations 
has a centered rectangular symmetry, c2mm. Its appearance is close to hexagonal, but the structure lacks the 7r/3 
rotational symmetry elements. If a hexagonal structure were formed from the wavevectors (0,1) and (\/3/2, 1/2), 
the ratio of the hexagonal spacing, Dh, to the lamellar spacing D would be 2/\/3 1.15. We find that the actual 
ratio of the lattice parameter of the stable, (non- aqueous), hexagonal phase to that of the metastable lamellar phase 
at the same temperature and composition to be 1.16. Thus the fluctuations do seem to be responsible for driving 
the lamellar phase to an ordered phase in which cylinders are oriented parallel to the lamellae, and with the correct 
spacing. The complete hexagonal symmetry, however, cannot be obtained from Gaussian fluctuations alone. We note 
that because the lamellae are isotropic in the x — y plane, all structure factors are isotropic in the k x — k y plane. 

For comparison, the same order parameter-order parameter structure factor, but now calculated at T* = 0.05, 
6 = 0.469, at which point the Gu phase is stable, is shown in Fig. 3. While the structure about the Bragg peaks 
remains, there are no satellite peaks here, and therefore no tendency to break the lamellae of headgroups into cylinders. 

The density-density structure factor, the Fourier transform of < <5$i(r)<5<i>i(r') >, with <5$i = —d<&2 the fluctuation 
in the total density of lipid 1, is obtained by setting the vectors a = b, defined in Eq. |3?| to be (1, 0, 0). As our two 
lipids are drawn from a homologous series and differ only in the length of the tails, this structure factor would not 
be easily measured. It is shown in Fig. 4 at T* = 0.0625, 8 = 0.356, and in Fig. 5 at T* = 0.05, 8 = 0.469. This 
structure factor is similar to the order parameter-order parameter structure factor. (One discernable difference is the 
value at zero wavevector which, for the density-density structure factor, is related to the osmotic compressibility. This 
small peak probably reflects the fact that it is rather easy to replace one lipid by the other because they have identical 
headgroups, and differ only by the length of their tails.) That the two different structure factors of Figs. 2 and 4 are 
so similar means that one of the two lipids tends to form cylinders, while the other must form their complement due 
to the constraint of incompressibility. Which lipid does which can not be determined from this correlation function, 
but it can from the cross correlation functions, as we now show. 

Consider the correlation between the density and order parameter fluctuations, < S^iS 1 ^ >. Because the headgroups 
of all lipids are identical, we expect that changes in the order parameter reflect, for the most part, changes in the 
densities of the tails. If the non lamellar-forming lipid 2 dominates the order parameter, then we expect that a 
positive fluctuation in the total tail density should be correlated with a positive change in the density of lipid 2. As 
the order parameter is defined as the density of heads minus the density of tails, a positive fluctuation in the tail 
density corresponds to a negative fluctuation in the order parameter. In addition, a positive fluctuation in in the 
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density of lipid 2 corresponds to a negative value of 5<f>i = —5<&2- As a consequence of the two minus signs, we expect 
< 5<S>i5ty > to be positive if the fluctuation in the non lamellar-forming tails dominate the order parameter. If they 
do not, it will be negative. A limit on this correlation can be obtained by noting that 



< (5<5> 1 ±Sy) 2 > >0. 



From this it follows that if 



then 



2 < <5$i<51< > 



< ((5$!) 2 > + < (6*) 2 >' 



1 > Ri > -1. 



(46) 



(47) 



(48) 



A value 1 of this ratio means that the density of lipid 2 is completely correlated with the order parameter, a value 
— 1 that the density of lipid 1 is completely correlated with the order parameter, and a zero value that there is no 
correlation between the density of either lipid and the order parameter. 

We have evaluated this quantity utilizing only the least stable fluctuation mode, i.e. that mode with non-zero 
wavevector whose energy will vanish at the spinodal. At T* = 0.0625, O = 0.356, near the lamellar spinodal, this 
wave vector is again k z — 1/2, k x ~ V3/2 . We find R\ ~ 0.27, which indicates that it is the nonlamellar forming 
lipid whose density is correlated with the order parameter, and that the magnitude of this correlation is appreciable. 

As a second measure of the importance of the nonlamellar forming lipid in the transition, we have considered the 
susceptibilities of each lipid to a disturbance with a wavevector equal to that of the least stable mode. At the spinodal, 
all responses at this wavevector will diverge, in general, but the amplitudes at which they do so will vary. Thus we 
have evaluated 



R 2 = 



_ < frfc 2 (k)<re 2 (-k) > 
< 5*i(k)<5*i(-k) > 



(49) 



at the wavevector of the least stable mode, where l (k) is the Fourier component of 6^ l (r) . At the same T* and 
9 near the lamellar spinodal, we find i? 2 = 2.05. Thus the susceptibility of the nonlamellar forming lipid to the 
perturbation at this wavelength is twice that of the lamellar forming lipid. 

A third measure of the relative importance of the two lipids is obtained by comparing the order parameter fluctua- 
tions of the two lipids to the order parameter in the Hu phase itself. We do this as follows. Consider the first star of 
wavevectors of the equilibrium hexagonal phase which is the stable one at the temperature and composition at which 
the lamellar phase is metastable. It consists of six wavevectors, kj, i = 1 to 6, of magnitude fc^) — Aiv/s/ZDh, where 
Dh is the lattice parameter of the Hu phase. The first four can be written as (k x , k z ) — fcC*)(±\/3/2, ±1/2), and the 
other two as (k x , k z ) — k^(0,±l). As noted earlier, if Dh/D were equal to 2/y/3, the first four wavevectors would 
lie on the edge of the Brillouin zone of the lamellar phase, while the others would correspond to reciprocal lattice 
vectors of the lamellar phase. We can determine the Fourier components with these six wavevectors of the order 
parameter fluctuations 5^2 and <5\&i . This tells us how large is the overlap of these fluctuations with the equilibrium 
hexagonal structure. We denote this Fourier component £\&i,(kj), L = 1,2. In general, it is complex so we calculate 
5*i(ki) + 5*i(-ki), which is real. At T* = 0.0625, 6 = 0.356, we find 



Rs 



^ 2 (k,) + ^ 2 (-k t ) 
<y*i(ki) + 5*i(-ki) 



1.43 
1.47 



i = 1 to 4 
i = 5,6. 



(50) 



If the fluctuations had hexagonal symmetry, the two values would be equal. We see that the amplitudes of the 
hexagonal fluctuations of the non lamellar-forming lipid are almost 50% larger than those of the lamellar-forming 
lipid. 

We next display the way in which the fluctuations of the lamellar phase near its spinodal indicate the beginning 
of the path to the stable Hu phase. We examine the effect of the fluctuation mode with lowest energy on the order 
parameter. In the absence of fluctuations, the order parameter, which is the difference in the total head and tail 
segment density, is 



J 1 ) A 2 ) 
11 ai n a 2 
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We add the contribution to this order parameter from the real part of the fluctuation mode with the lowest energy, 
i?e{(54'(r)}. We add it with a variable amplitude, e to better visualize its effectEil. Fig. 6(a) shows the order parameter 
in mean-field theory, with no fluctuation contribution, e = 0. The point in the phase diagram is T* — 0.0625, 9 = 0.356. 
The center of the headgroup region is at z — 0. The lamellae are in the x — y plane, and the coordinates are measured 
in units of the lattice spacing, D. The grey scale divides the values from —0.6 to 0.6 into ten shades. The darkest 
values are most positive, and correspond to a dominance of head groups. Values in (a) range from —0.214 to 0.251. 
In (b), (c) , and (d), the fluctuation contribution is turned on with amplitudes e = 0.06, 0.12, and 0.18 respectively. 
The extreme values in (d) are —0.301 and 0.572. One clearly sees the undulation of the tail region, and the pinching 
off of the head group region, leading to an Hjj phase in which the location of the tubes is coplanar with the lamellar 
phase. Again, the lattice spacing of the stable Hu is within 1% of the spacing expected if there were a perfect epitaxy 
between the phases. The mechaaism we see indicated here by the fluctuations from the lamellar phase is very similar 
to those proposed by Hui et al.ll3 and Caffreyia. The coplanarity which is a direct result of it has been observed by 
Gruner et al. with orientecL-nhotoreceptor membranesEj. We also note the possible significance for this problem of 
an observation made in RefBi. Due to isotropy in the x — y plane, the fluctuation modes are infinitely degenerate in 
the kx 

plane, that is, all fluctuation modes with the same magnitude k\ + ky have the same energy of excitation. 
If one particular direction is preferentially excited, the mode leads to the formation of ordered cylinders as we have 
seen above. But if several directions are excited simulataneously, then the ripples in the corresponding directions will 
interfere with one another. .Qn<? pattern that can be obtained is a periodic array of holes,, in the lamellae, pmdncing 
a perforated lameller phascE3E3. These perforations are similar to pores in the membraneE9Ed, and to stalksQOEl. 

Finally, we display the effect of the fluctuation mode with the lowest energy on the difference in the tail distribution 
of the two lipids. The tail distribution of lipid L as calculated in the self-consistent field theory, without fluctuations, is 
0t z,(r). We add the contribution to this distribution from the fluctuation mode with the lowest energy, eRe{5$> t ,L( r )} 
where e is, again, an amplitude which we can adjust. In order to compare the tail distributions of the two lipids, 
we must take account of the fact that, because the concentrations of the two lipids are unequal, they have different 
average tail segment densities 



J t,L 



1 

V 



dv &, L (r). (51) 



Therefore we utilize the quantity 



and plot 



P tL (e,r)^^ F)+ f e{S ^ } (52) 

4>t,L 



AP t (e, r) = P tll (e, r) - P t , 2 (e, r) (53) 

for various values of e. This shows us how the fluctuations about the lamellar phase tend to spatially arrange the 
two different lipids. Again we place ourselves at the point in the phase diagram at which T* = 0.0625, = 0.356. 
Figure 7(a) shows this difference within the self-consistent field theory alone, that is with e = 0. The center of the 
headgroup region is at z — 0. The grey scale is such that darkest values correspond to positive values of Pt at which 
lipid 1 dominates, and lightest values correspond to negative values of Pt at which the non lamellar- forming lipid 2 
dominates. Thus the density of lipid 2 is largest at the center of the tail region, which is expected as lipid 2 has the 
longer tails. The maximum and minimum values on the grey scale of Fig. 7 correspond to —0.16 and 0.16, divided 
into ten shades of grey. The maximum and minimum values attained by AP t (0,r) are 0.089 and —0.109. Thus the 
differences in local tail segment densities shown in Fig. 6 are on the order of 10% of the average tail segment densities. 

In figs. 7(b), (c), and (d) we turn on the contribution to the density difference of the tail segments from the 
lowest energy fluctuation mode, and plot APt(0.06,r), APt(0.12,r), and APt(0.18,r). The largest and smallest 
values attained by AP t (0.18,r) are 0.160 and —0.147. In this sequence, one clearly sees that in the fluctuations, the 
non lamellar-forming lipid dominates the undulations of the tail region, particularly in the direction of the second 
neighbors, the direction which is most difficult for the tails to fill. The tails of the lamellar-forming lipid are initially 
relegated to filling in the region between the cylinders which form within a given lamellae. 

In sum, we have employed a model system of a non-aqueous mixture of lamellar- and nonlamellar-forming lipids to 
examine the Gaussian fluctuations of the lamellar phase near its spinodal. We have calculated the structure factor 
which would be experimentally measurable were the lameller phase sufficiently oriented. We have found that the initial 
stage of the L a to Hu transition involves the formation of undulations in the tail-rich lamellae, undulations which 
are dominated by the nonlamellar-forming lipid. As we have noted, the lattice constant of the ordered, stable, Hu 
phase is within 1% of that which would be produced by the lowest energy fluctuation mode, and the orientations are, 







of course, identical. There is also a coplanarity between one of the principal direction of the tubes in the H a and the 
orientation of the L a phase which produced it, a coplanarity observed in experiments. We have used several measures 
to show that the nonlamellar-forming lipid plays a dominant role in this transition in the mixture. Because we have 
examined the fluctuations very near the spinodal, we are presumably seeing the beginning of spinodal decomposition. 
Further from the spinodal, and in particular, near the L a phase boundary which is generally far from the spinodal 
(see Fig. 1), it is to be expected that the transition proceeds by a nucleation and growth mechanism, which may well 
involve some form of lipid intermediate structureoc3. Presumably by making a series of controlled quenches from 
the lamellar to the hexagonal phase, one should observe the progression from the spinodal decomposition mechanism 
of undulations to the nucleation and growth intermediates. Certainly both mechanisms seem to have been observed 
in the same systemLJ. In either form of the transition, structures which are nonlamellar must be created. As we have 
shown and made quantitative, the nonlamellar forming lipids tend to dominate the process of this kind which we have 
examined. This lends support to the line of argument that their role in the biological membrane is to facilitate the 
creation of such structures. 

This work was supported in part by the National Science Foundation under grant number DMR9876864. 
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Figure 1 Phase diagram of an anhydrous mixture of lipids 1 and 2 in the plane of temperature, T* , and volume 
fraction, 9, of lipid 1, the lamellar- forming lipid. Phases shown are inverted hexagonal, Hjj, inverted gyroid, 
Gn, lamellar, L a , and disordered, D. The dashed line is the spinodal of the L a phase. 

Figure 2 Order-parameter order-parameter structure factor shown in the k z — k x plane. The wavevectors are in units 
of 2ir/D, with D the lamellar spacing. The point in the phase diagram, T* — 0.0625, = 0.356, at which the 
factor is calculated, is near the spinodal of the lamellar phase. 

Figure 3 Same structure factor as in Fig. 2, but at the point in the phase diagram T* = 0.05, = 0.469, far from 
the spinodal of the lamellar phase. 

Figure 4 Density density structure factor calculated at T* — 0.0625, = 0.356. 

Figure 5 Density density structure factor calculated at T* — 0.05, = 0.469. 

Figure 6 The total order parameter is shown in the x— z plane with different amplitudes, e, of the contribution from 
the fluctuation mode with the lowest energy; (a) e = 0., (b) e = 0.06, (c) e = 0.12, (d) e = 0.18. Light areas 
indicate regions in which tail density is greatest, dark regions where head density is largest. 

Figure 7 The difference in the tail segment densities arising from lipid 1 and from lipid 2. The amplitude of the 
contribution from the fluctuation mode with the lowest energy is (a) e = 0., (b) e = 0.06, (c) e = 0.12, (d) 
e = 0.18. Light areas indicate regions in which tails of lipid 2, the nonlamellar former, dominate, dark regions 
the predominance of segments from the tails of lipid 1. 
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FIG. 2. X.-j. Li, J. Chem. Phys. 
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FIG. 5. X.-i. Li. J. Chem. Phvs. 



